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We explore theoretically the magnetoresistivity of three-dimensional Weyl and Dirac semimetals 
in transversal magnetic fields within two alternative models of disorder: (i) short-range impurities 
and (ii) charged (Coulomb) impurities. Impurity scattering is treated using the self-consistent Born 
approximation. We find that an unusual broadening of Landau levels leads to a variety of regimes of 
the resistivity scaling in the temperature-magnetic field plane. In particular, the magnetoresistance 
is non-monotonous for the white-noise disorder model. For H ^ 0 the magnetoresistance for short- 
range impurities vanishes in a non-analytic way as . In the limits of strongest magnetic fields 
H, the magnetoresistivity vanishes as l/H for pointlike impurities, while it is linear and positive in 
the model with Coulomb impurities. 


I. INTRODUCTION 

Topological materials and structures represent one of 
the central research directions in the modern condensed 
matter physics. One of the classes of such materials is 
topological insulators and superconductors which possess 
a bulk gap and topologically protected surface excitations 
with massless Dirac spectra. Another important class 
is gapless materials with topologically protected Fermi 
points or nodal lines. The most well-known example is 
graphene whose dispersion is characterized by two Fermi 
points where the valence and the conduction bands touch. 
Excitations in the vicinity of these points have a linear 
dispersion and can be viewed as two-dimensional Dirac 
fermions. 

Three-dimensional counterparts of graphene are Dirac 
semimetals. In such a material the dispersion near the 
nodal point is characterized by 4 x 4 Dirac Hamiltonian, 
i.e. the conductance and the valence bands have an ad¬ 
ditional twofold degeneracy. Experimental realizations 
of Dirac semimetals include Cd 3 As 2 ^ and NasBi^; fur¬ 
ther candidate materials have been recently discussed^. 
The twofold degeneracy discussed above can be lifted if 
either spatial inversion or time-reversal symmetry is bro¬ 
ken. The four-component Dirac solution then decouples 
into two independent two-component solutions represent¬ 
ing two Weyl fermions of opposite chirality. Thus, each of 
the Dirac points then splits into two Weyl points without 
any additional degeneracies. Recent experiments pro¬ 
vided evidence that TaAs^’® and NbAs® can be classified 
as Weyl semimetals^’®. Further promising candidates for 
Weyl semimetals include pyrochlore iridates®, topologi¬ 
cal insulator heterostructures^°, and Cd 3 As 2 with low¬ 
ered symmetry^^. In the rest of the paper we will use the 
term “Weyl semimetal” in a broader sense, including also 
the degenerate case of Dirac semimetals. 

Transport properties of Weyl semimetals are highly pe¬ 
culiar (for various theoretical aspects of the problem, see, 
e.g.. Refs. 9, 12-27 and references therein). In particu¬ 
lar, weak disorder (with a strength below a certain criti¬ 
cal value) has a negligible effect on the density of states. 


Specifically, the density of states vanishes quadratically 
in energy around the Weyl point despite the presence 
of disorder. Furthermore, the limits T —> 0 and w —)■ 0 
(where uj is the frequency) are not interchangeable for the 
behavior of the conductivity (assuming that the chemical 
potential is at the Weyl point). While sending frequency 
to zero first results in a finite conductivity, the zero-T 
ac conductivity vanishes®’^^’^^’^® in the limit a; —>■ 0 as 
|w|. In the strong-disorder regime these singularities are 
eliminated. 

Behavior of Weyl semimetals in the external magnetic 
field is also expected to be very nontrivial. This is re¬ 
lated, first of all, to the unconventional Landau quan¬ 
tization of Dirac fermions. Furthermore, since a sin¬ 
gle Weyl node displays a chiral anomaly, a possibil¬ 
ity to control the valley polarization as well as a large 
anomalous Hall effect are expected.Much atten¬ 
tion has been recently put on the longitudinal magne¬ 
toresistance in Weyl semimetals that originates from the 
chiral anomaly. 

In this paper, we develop a theory of the transversal 
magnetoresistivity of a Weyl semimetal. We are par¬ 
ticularly interested in the range of sufficiently strong 
magnetic fields, such that the Landau quantization is 
important.®® One of the motivations for our work was a 
spectacular experimental observation of a large, approx¬ 
imately linear magnetoresistance in the Dirac semimet¬ 
als Cd 3 As 2 and TlBiSSe in strong magnetic field. 
Quantum linear magnetoresistance has been obtained by 
Abrikosov in a seminal paper. Ref. 39, for Dirac semimet¬ 
als in the extreme limit when only one Landau level 
is filled (i.e., the cyclotron frequency exceeds the tem¬ 
perature). The linear behavior was traced back to the 
magnetic-field dependent screening of charged (Coulomb) 
impurities. 

We consider a general case of arbitrary relation be¬ 
tween the magnetic field and temperature such that, de¬ 
pending on the regime, the magnetoresistance is dom¬ 
inated by contributions from the zeroth Landau level, 
other separated Landau levels, and overlapping Landau 
levels. We address two alternative models of disorder: (i) 
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short-range impurities and (ii) Coulomb impurities. We 
show that for pointlike impurities the transition between 
the weak-disorder and the strong-disorder phases persists 
in the presence of magnetic field and explore singulari¬ 
ties of the weak-disorder phase. Further, we find that an 
unusual broadening of Landau levels leads to a variety 
of regimes of the resistivity scaling in the temperature- 
magnetic field plane. The transversal magnetoresitance 
is found to be non-monotonous for the model of weak 
white-noise disorder. Remarkably, in the limit of 7L —>■ 0, 
the magnetoresistance for short-range impurities shows a 
non-analytical behavior For strong pointlike im¬ 

purities and for charged impurities we find a positive 
quadratic magnetoresistance at iL —>■ 0. In the limit of 
strongest magnetic fields H, the magnetoresistivity van¬ 
ishes as 1/iL for pointlike impurities, while it is linear and 
positive in the model with Coulomb impurities, in agree¬ 
ment with Abrikosov’s result^® and recent experimental 
findings. 

The paper is organized as follows. In Sec. II, we out¬ 
line the implementation of the Born approximation in 
the context of Weyl semimetals subjected to a magnetic 
field in the presence of pointlike impurities. We find that, 
in analogy with the zero-magnetic-field case, one has to 
distinguish between weak and strong disorder regimes 
separated by a phase transition. In Sec. Ill we develop 
the formalism of the self-consistent Born approximation 
(SCBA) and analyze the broadening of Landau levels due 
to disorder and magnetic field. Section IV presents the 
formalism for calculation of the conductivity in magnetic 
field for the model of white-noise disorder. In Sec. V, we 
extend our analysis to the case of charged impurities. In 
Sec. VI we use the obtained results to calculate and an¬ 
alyze the magnetoresistance for both models of disorder. 
Our findings are summarized in Sec. VII. 


II. LANDAU-LEVELS SPECTRUM IN WEYL 
SEMIMETALS 

A. Clean case 

In the presence of a constant homogeneous magnetic 
field H in z direction, electrons of a single Weyl node are 
described by the Weyl Hamiltonian 

(p) = J d^r’^\r)v(T (^p --A^ (1) 

where cr denote the Pauli matrices and A(r) = (0, Hx^ 0) 
is the vector potential. (We have chosen the Landau 
gauge.) The eigenfunctions of this Hamiltonian have 
two components (a,/3 = 1,2) in the space spanned by 
(T. Positions of the Landau levels (LLs) in a clean Weyl 
semimetal are given by 


£0 =vpz, 

(2) 

= ±vJpl + 

(3) 



FIG. 1. Density of states of a clean Weyl semimetal. The blue 
solid line shows the DoS in a finite magnetic field, Eq. (7), 
while the red dashed line corresponds to the zero magnetic 
field, Eq. (6). 


where Ih = {sH/c) is the magnetic length and 
n = v^^eHjc is the distance between the zeroth and 
first LL. We set h = 1 throughout the paper. In the fol¬ 
lowing, we choose the energy bands such^® that the wave 
function of the clean zeroth LL has only component 1 in 
the pseudospin space: 


= 


72 



1/2 gj(PH!/+PiZ) 
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1/2 ^iiPyV+Pzz) 

z 


1 ( 3 : 


l^Py) 

(4) 


for n > 0 and = 9{±pz)<^o, = 0 for n = 0. 

Here (pn are the normalized eigenfunctions of free elec¬ 
trons in magnetic field and 6 denotes the Heaviside step 
function. The retarded bare Green function Go of the 
clean system is conveniently represented as a matrix in 
the pseudospin space of bands a, f3 = 1, 2: 


G^°l= 


'i’na 'i’n/3 


^ £ -I- to - e7 

n>0,X—± ^ 


(5) 


It is worth noting that the summation of A eliminates the 
theta-functions 6{zLpz) in the n = 0 term in Eq. (5), so 
that the integration over p^ in what follows will always 
be performed from —00 to 00 . 

In the zero magnetic field the clean density of states 
(DoS) behaves quadratically in energy: 


u{e) 




( 6 ) 


In the presence of magnetic field, the DoS acquires a 
sawtooth form with square root singularities originating 
from the one-dimensional {pz) dispersion of each Landau 
band (cf. Ref. 40) 


/(£) = 




1 + 2 


e+l/2 
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7 ^^ - 2nv^/ll 
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(7) 
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The DoS is visualized in Fig. 1. 


B. Introducing disorder 

We consider now the effect of disorder. The im¬ 
purity scattering generates a self-energy E(p,£) in the 
(impurity-averaged) Green function, 

^ e - uo- • (p - £A) - S(p,e)’ 

, , ( 8 ) 

which is a matrix in the pseudospin space (in which the 
Pauli matrices cr operate). 

We will assume that the disorder potential is diagonal 
in both spin and pseudospin indices and neglect scatter¬ 
ing between different Weyl nodes. We will discuss this 
approximation and its limitations in the end of the paper, 
Sec. VII. Clearly, in the absence of internode scatter¬ 
ing, the structure in the node space will be trivial for all 
quantities. We do not show it explicitly below; the cal¬ 
culated density of states and the conductivities are those 
per Weyl node. 

We will first consider a model of pointlike impurities 
and later analyze generalization to the case of Coulomb 
impurities. The impurity potential has then the form 

Vdis(r) = ■uo^(5(r - Ti)!, (9) 

i 


where 1 is the unit matrix in the pseudospin space. In 
view of the matrix structure of the impurity potential 
Vdis(r), the impurity correlator W becomes a rank-four 
tensor. The self-energy reads 



d^q 

(27r)3 




( 10 ) 


read: 

r = _ £ - S 2 + VP^ _ . . 

(£ - Si - UPz)(£ - E 2 + WPz) - 

^ _ _ £ - Si - vp^ _ 

(e - Si - vpz){e - S 2 -I- vpz) - Vt‘^{n + 1)' 

(14) 

In general, the self-energy depends on energy and on the 
LL index, S = T,(e,pz,n). However, for a white-noise 
disorder, the dependences on n and p^ drop out. 


C. Born approximation 


We start with the Born approximation, where we ne¬ 
glect the self-energies in Green’s functions (13) and (14) 
for the calculation of self-energies: 


Sf(£) 


S?(e) 


7 f°° <ip^ _ £ + vPz _ 

27rl|i ^ Stt (£ -F i0)2 - H2n - v'^pl ’ 

7 ^ _ £ - vPz _ 

^ 27r {e + iO)"^ - + 1 ) - v'^pl' 

(16) 


The summation over n here should, in fact, be restricted 
by an upper cut-off iVinax, as will be discussed below. 
After shifting the summation over n in S 2 , we see that 
the two self-energies differ only by the absence of the 
n = 0 term in S 2 : 


7 dpz £ + vpz _ 

27rZ|j 7-00 27r {e + lO)^ - v^p^ ~ 


iA + 


where, using I'jj = 2 (u/H)^, we have introduced 


2 A£ 
ttA ’ 
(17) 


Sttv^ ’ 


(18) 


For a diagonal impurity potential, the impurity correlator 
is diagonal as well, which is expressed as 


kFQ,^^,5(q) — 7^07^/3(55 (^^) 

where 7 = riimpUg. The self-energy is diagonal in the 
energy-band space. However, in the presence of magnetic 
field, the self-energy is no longer proportional to the unit 
matrix: 


S = diag(Ei,E 2 ). ( 12 ) 

This asymmetry originates from the asymmetry of states 
in the zeroth LL. In the clean case, the states of the 
zeroth LL are only present in one energy band. Later we 
will see that a strong impurity scattering eliminates this 
asymmetry. 

We switch to LL representation so that G = G(e,pz, n) 
and E = 'E(e,pz,n). The diagonal components of the 
matrix Green function ( 8 ) that determine the self-energy 


and A is the bandwidth. The necessity to introduce the 
ultraviolet cut-off A originates from the approximation of 
a true energy dispersion by the Dirac-fermion one, which 
is, in fact, a low-energy approximation. Our analysis is 
applicable for £, S(£) <C A. 

In view of Eq. (17), it is sufficient to calculate Si: 

ImE.(r) = -A|£| ( 19 ) 

n—0 

Njnax 

ReSi(£) = —Ae ^ 

n—A^e + 1 

The imaginary part of the Born self-energy is produced 
by the Landau levels below s, while the real part is due 
to the contribution of the Landau levels above £. In 
Eqs. (19) and (20) 


— £2 


( 20 ) 


Ne = 



( 21 ) 
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is the number of the Landau level below energy e, and the 
symbol [...] denotes the integer part of a number. For the 
evaluation of the sum in Eq. (20), we have introduced the 
upper cutoff iVmax which is determined by the ultraviolet 
energy cutoff A in the following way: A^max = 

(A^max is the index of the highest Landau level within 
the bandwidth A). 

The sum in ReSi is dominated by the upper limit 
A^max^ ReEi ~ e(A/n)Afm4x- Assuming no Landau 
quantization at the ultraviolet energies ^ A, we use the 
zero-iL result 


ReEi(£) 



with 


P 


7 A 


( 22 ) 


(23) 


The parameter P quantifies the strength of disorder. For 
sufficiently strong disorder, the real part of Born self¬ 
energy becomes larger than e, which clearly signifies the 
insufficiency of the simple Born approximation. As we 
will see in Sec. IIIC below, the self-consistent treatment 
of strong disorder yields a dramatic change of the behav¬ 
ior of the density of states for strong disorder. 

In what follows, however, we mostly focus on the limit 
of weak disorder, /3 <C 1. We absorb ReE into the redefi¬ 
nition of the energy e ^ e = £{1 + P/2), and neglect the 
difference between e and e. For |e| < fl, we find 


ImSi = -A, ImEa = 0. (24) 


For higher energies, using the Euler-Maclaurin formula 
for the sum over n < — 1, we express the imaginary 

part of the Born self-energy as 


ImEi(£) ~ —A 
1 


Ve^ - V ye) 


1 + 


v /£2 _ (N, - l)fI2 


2 ^ 


(25) 


This result is illustrated in Fig. 2. The first term in 
the square brackets of Eq. (25) is responsible for the 
square-root divergency at the positions of LLs, whereas 
the last term yields the parabolic background similarly 
to the zero-iL case. 


III. SELF-CONSISTENT BORN 
APPROXIMATION 

We now turn to the self-consistent Born approxi¬ 
mation (SCBA). The self-consistent treatment is moti¬ 
vated by the presence of square-root singularities in the 
Landau-level broadening (25) obtained within the Born 
approximation. The introduction of disorder-induced 
self-energies in Green’s functions should cut off such di¬ 
vergencies. 


ImZ] 


A 



FIG. 2. Disordered Weyl semimetal in Born approximation. 
ImEi(£) in units of A. Blue: Eq. (19); Dark blue, dashed: 
Eq. (25); Green: 2{£/Q2p — ie/2Q,-\-l/2 [minima in Eq. (25)]; 
Red: 2(e/f2)2 (result without magnetic field). 


The SCBA equations (10), (13) and (14) with the dis¬ 
order correlator (11) acquires the form (below z = vpz). 


Ei(£)=A^ / dz 


€ — S 2 z 


n>0 ^ 


(e — Si — z){s — S 2 + z) — 

(26) 


E2(£) = A ^2 / dz 


£ — El — z 


n>l ‘ 


(£ — El — z)(£ — E 2 -h z) — 

(27) 


As above, we absorb the real parts of self-energies (de¬ 
termined by the ultraviolet cut-off A) into the shifts of 
energies £—>■£. The density of states is given by 


p(e) = -iTr ImG = —- (ImEi -h ImE 2 ). 

TT 7r7 


(28) 


In the regime of well separated LLs (the corresponding 
conditions will be analyzed below), the sum over the Lan¬ 
dau levels is evaluated as follows. Let us assume that the 
energy is close to the bottom of the N-th Landau level 
Ng ~ N. Then the main contribution of the sum over 
n comes from the term n = N. We thus single out the 
term oi n = N from the sum and evaluate it separately 
from the sum over the remaining Landau levels. 

We note that the self-consistent treatment of the LL 
broadening is fully justified for weak disorder and e 
fl. All the renormalization effects not captured by the 
SCBA affect the real part of the self-energy, see discus¬ 
sion in Ref. 41 where the SCBA was employed for 2D 
Dirac fermions in graphene. In that case, disorder was 
marginally relevant and its effect could be incorporated 
through the renormalization of parameters (induced by 
the contributions of higher LLs) that enter the SCBA 
equations for a given Landau level. In the present 3D 
case, the renormalization of ReE can be neglected for 
the case of weak disorder, /3 1. Even for the lowest 

LL, the SCBA density of states in two dimensions is para¬ 
metrically correct.The exact shape can be calculated 
along the lines of Refs. 42 and 43 In three dimensional 
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systems for weak disorder, the extra integration over the 
momentum pz further reduces the difference between the 
exact and SCBA results in the limit £ —>■ 0. 


A. Energies close to the lowest Landau level 


We first consider the case of lowest energies, |£| <C 
(i.e., W = 0) for weak disorder, /3 <C 1. In this case, the 
asymmetry with respect to the zeroth LL should be taken 
into account and the imaginary parts of the two self¬ 
energies strongly differ from each other. When the lowest 
Landau level is well separated from the others (ImSi _2 'C 
n ), the contribution of higher Landau levels to the sum 
over n can be treated within the Born approximation, 
while the contribution of n = 0 should be calculated self- 
consistently. Within this procedure, we immediately get 
Iml ]2 = 0 and 


/ OO 

dz 

-OO 




Iml]i(£) 

zY + [Iml]i (£)]2 


-A. 

(29) 


This result coincides with the result of the non-self- 
consistent Born approximation, Sec. IIC. The zeroth 
LL is separated from the first as long as the condition 
A < n is fulfilled. The density of states for £ ^ 17 is 
finite and, to the leading order, is energy-independent. 

Using Eq. (29), we find the leading non-vanishing term 
in Iml ]2 for £ ^ 17: 


ImE 2 (£) 



n—1 


fi^n 

(fi^n — £2^3/2 


-A/3. (30) 


Thus, in the limit of weak disorder, /3 1, we can neglect 

ImE 2 for £ ^ 17. In fact, the condition is even softer: 
ImE 2 becomes of the order of ImEi only in the close 
vicinity of the first Landau level, |£ — 17| ^ A. 


B. Energies at high Landau levels 


We now consider high energies, £ ^ 17. As we have 
seen in Sec. IIC, already within the Born approxima¬ 
tion the average (as well as minimal) broadening of Lan¬ 
dau levels increases with e parabolically, as in the zeio-H 
case: oc A{e/flY ^ see Fig. 2. Therefore, the 

difference between ImEi and ImE 2 that comes from the 
contribution of n = 0 can be neglected for energies e away 
from the zeroth LL. In what follows we set Ei = E 2 = E 
for £ ^ 17. 

Introducing 


ri,2(£) = -ImEi,2(£) (31) 

and setting Ti = r 2 = T, we arrive at the self-consistent 


r 

A 



FIG. 3. Disordered Weyl semimetal in the self-consistent Born 
approximation for e ^ 17: r(£) in units of A obtained by 
numerical solution of Eq. (32). Blue, red, and green curves 
correspond to A/17 = 10“'*, 10“^, 10“^, respectively. For all 
curves Amax = 100. 


equation for r(£ ^ 17) : 


r = ^r(")(£), 

n—0 


(32) 


_ at £2 -b n^n + r2 

TT 7_oo ^ (£^ — 17^n — — 2:^)2 -b 4 e 2 p 2 


= A Re 


f£ -b r 


(33) 


where we have introduced the partial contribution T*^”) (£) 
of the nth Landau level to the total broadening r(£). 
Note that each term in the r.h.s. of Eq. (32) contains the 
total broadening T rather than the partial T*^”^. Further, 
the position of the nth Landau level is shifted by disorder: 
17^n appears only in combination 


W^ = n^n + T^. (34) 


In the case of weak disorder, for all energies 17 <C £ <C A 
we have £ 7^ r(£), so that Eq. (33) can be written as 


r(”)(£) ~ A£ 


^£2 - + ^( W 2 - £ 2)2 + 4 £ 2 r 2 

y/(IU 2 _ 5 : 2)2 + 4£2r2 ' 


(35) 


For e ^ OO this yields r*^”^(£) A and for e Wn we 
get r("^(£) —>• AT/Wn- Thus, when £ crosses Wn, the 
nth Landau level gets an extra contribution A to r(£). 
The solution of the self-consistent equation (32) is shown 
for A/17 = 10“4, 10“3, 10“2 in Fig. 3. 

Let us now fix the Landau-level number N ^ 1 and 
consider the range of energies around Wn- Assuming 
well separated Landau levels below N, we neglect T in 
all terms with n < N: 


N-l 


N-1 


^r(")(£)c£A^ 


n—0 


pi _ W 2 
n=0 n 


~ 2AA. (36) 
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The contribution of Landau levels with n > fV+1 is dom¬ 
inated by A^max and can be neglected for weak disorder: 


N„ 




WdT 




n—N-\-l 


(W^ - £2)3/2 


r/ 3 <r. 


(37) 

Finally, the contribution of the fVth Landau level (closest 
to the energy e) can be further simplified for \£—Wm\ 'C 
Wn-. 


r(^)(£) ~ 


^s-WN + ^/{WN-e)^ + T^ 


- eY + r2 


In particular, exactly at £ = Wn we find 

^£ 1/2 


rW(£ = VF^)~ 


2 ri /2 • 


(38) 


(39) 


Using Eqs. (36) and (38), when £ is close to Wn, the 
self-consistency equation takes the form: 


2 £l £2 

^ ^ ^ ^ - £)2 + r27^ ’ 

_ (40) 

where rCg ~ ily/N^, so that the r.h.s. of Eq. (40) explic¬ 
itly depends on e only, as it should be. 

Exactly at £ = Wn, the self-consistency equation 
reads: 


2yl£2 yl£^/2 

U2 2ri/2 ■ 


(41) 


We observe that for sufficiently small energies, the broad¬ 
ening is dominated by the self-consistent contribution of 
the same Landau level, whereas for large energies, the 
broadening is given by the zeio-H result stemming from 
lower Landau levels: 


r(£ = Wn) 


(A/2)2/3el/3^ U<£<£*, 

2 A(£/U) 2 , £>£*, 


where 


112/5 

£* ~ n{n/A)^/^ (X 

/yi /0 


(43) 


Below £* Landau levels are fully separated. Each peak in 
r(£) is non-symmetric with respect to Wn, as inherited 
from the clean density of states. The shape of the LL 
broadening is analyzed in detail in Appendix A. 

For high energies, the behavior in zero magnetic field 
should be recovered. Indeed, we can express the result 
for £ > £* in terms of the energy as follows: 


r(£) = ^£ 2 . ( 44 ) 

47ru'^ 


into account the corrections to the broadening at £ > £*, 
we will see in Sec. HID below that the Landau level 
quantization of the density of states remains intact in 
a finite range of energies above e*. This should be con¬ 
trasted with the 2D case, where a single scale separates 
regimes of strong and weak Landau quantization. Fi¬ 
nally, we note that magnetooscillations in Weyl semimet¬ 
als were addressed in Ref. 44 with phenomenological 
energy-independent broadening. We find, however, that 
the energy dependence of F is very rich. 


C. Strong Disorder 


We now briefly discuss the regime of strong disorder, 
/? > 1. As follows from the consideration of the weak- 
disorder case, see Eq. (30), for strong disorder the dif¬ 
ference between the two self-energies. Si and E 2 , be¬ 
comes inessential even at £ = 0. After the evalua¬ 
tion of the sum over n in Eqs. (26) and (27), we find 
a qualitative change in the behavior of the imaginary 
part of self-energy (and thus of the density of states) at 
AAy/Nraax = This implies the existence of a criti¬ 
cal disorder strength, jc,mag = 2 tt^v^/A separating the 
two regimes. In the absence of magnetic field, the emer¬ 
gence of such a critical disorder strength 7 c was reported 
in Refs. 16, 17, 45, 21 and 24. Remarkably, the criti¬ 
cal disorder strength ^c,mag which we find for the case 
of a strong magnetic field turns out to be equal to the 
zero-field value 7 c. 

The solution of the SCBA equation for strong disorder 
for £ ^ U is given by 


r ~ 


2Uv'f/Vmax 

TT 


02 

2 A 




(45) 


which is equal to the zero-H result obtained in Ref. 16. 
When 7 is substantially larger than 7 c (i.e., 7 —7c ^ 7c), 
the broadening becomes of the order of the ultraviolet 
cut-off. 


F - A, (46) 

which ensures that all Landau levels overlap. Further, at 
£ <C U, the solution of the SCBA equations yields for the 
LL broadening 

ImSi ~ ImE 2 ~ -Oa/ N^ax ~ A. (47) 

Thus, when the disorder is substantially stronger than 
the critical one, even the zeroth LL overlaps with the 
rest of the spectrum. 

Within the SCBA, the real part of the self energy for 
P = ihc 1 is found to be 

ReSi ~ ReE2 ~ |^£, (48) 


We see that the magnetic field has dropped out from this 
result, as expected. Thus the LL broadening is domi¬ 
nated by the H = 0 result for £ < £* In fact, taking 


yielding £ ~ e/p. 

In this paper, we do not discuss the critical regime 
near the transition from weak to strong disorder at /3 ~ 
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FIG. 4. Density of states, vie), within self-consistent Born 
approximation in units of v{Q) oc H [Eq. (52)], as obtained 
from Eq. (32) for e > and Eq. (29) for e < 0,. The curves 
corresponds to (a) A/Q. = (b) A/Q. = 10“^, and (c) 

A/^ = 10“^. The value A^max = 100 was used. 


1. At zero magnetic field, the criticality was addressed 
using the e-expansion within the renormalization group 
approach in Refs. 20 and 24. The effect of magnetic field 
near the transition remains a very interesting question 
for future work. 


D. Density of States 

In Fig. 4, we plot the density of states obtained by a 
numerical solution of the SCBA equation in the case of 
weak disorder. The three figures illustrate the evolution 
of the density of states with the increasing value of the 
parameter (proportional to the disorder strength 

and to the square root of the magnetic field). 

Even in the clean case, Landau levels are broadened 
due to the integration over see Fig. 1, so that the 


divergent peaks at £n{Pz) are located on top of the back¬ 
ground density of states. However, these peaks are well 
resolved for all energies. Disorder leads to the suppres¬ 
sion of the peaks and eventually Landau levels fully over¬ 
lap at high energies. Let us discuss the characteristic val¬ 
ues of Landau level index at which the behavior of the 
density of states changes qualitatively. 

At £ > £*, the broadening of Landau levels is domi¬ 
nated by the background (zero-77) contribution. From 
Eq. (43), we see that the corresponding LL index AL = 
decreases with increasing H: 

N. ~ (n/Af/‘ « (49) 

In order to find out whether LLs are resolved or 
not, we should check whether the condition Sn+iiPz = 
0) — EniPz = 0) > F is fulfilled. For energies D ^ £ <C 
n(n/A)^/®, the width of high Landau levels is smaller 
than the distance between them. Therefore, at 

£ ~ £* we have a situation when Landau levels are still 
resolved on top of the background, but the height of the 
small peaks , Eq. (41), is lower than the height 

of the background. The latter then dominates the broad¬ 
ening. 

The neighboring peaks fully overlap in the regime £ > 
£** > £*, where the broadening is given by the zero-77 
result: 

£2 7)2 

V-T ^ = 

The corresponding LL index 

~ W ~ \ a) 72/3pfl/3 ) 


decreases with increasing magnetic field, similarly to N^,. 
Let us emphasize that, contrary to conventional expec¬ 
tations, the number of separated LLs decreases with in¬ 
creasing magnetic field. We thus see that a large number 
of low-lying Landau levels are well resolved for weak mag¬ 
netic field. 

When the magnetic field increases for a fixed disorder 
strength, the Landau level index AL* associated with the 
starting point of overlapping, becomes smaller. This be¬ 
havior is very unusual, as it is opposite to that in the 
case of conventional semiconductors. The energy £** 
where the Landau levels start to overlap increases with 
77 as 77^/^. With regard to the behavior with disorder 
strength, the obtained results qualitatively conform with 
intuitive expectations. Specifically, with increasing dis¬ 
order, the number of separated Landau levels decreases 
and the corresponding energy range shrinks. 

It is important to stress that, in the presence of mag¬ 
netic field, the density of states at the Dirac point (zero 
energy) is finite even in the weak-disorder regime. Specif¬ 
ically, the value of the density of states at £ = 0 is linear 
in magnetic field: 


^( 0 ) 


A 

Try 


7)2 


(X 77. 


(52) 





























separated LL’s, 
but larger 

zeroth LL separated LL’s background overlapping LL’s 



n 

(7fi)l/5 


FIG. 5. Relevant energy scales and dominant contributions to the density of states of a disordered Weyl semimetal (for brevity, 
we have set v = 1). 


It is worth noting that a finite value of the density of 
states at the degeneracy point will lead to a finite con¬ 
ductivity independently of the order of limits w —>■ 0 and 
T —)■ 0. From this point of view, a finite magnetic field 
has the same effect as a strong disorder. 

All the features of the density of states that we have 
found analytically (see Fig. 5) are perfectly observed in 
Figs. 4 a,b,c. First, one sees that for weak disorder and 
weak magnetic field many LLs are separated and that 
the number of separated Landau levels decreases with in¬ 
creasing magnetic field or with increasing disorder. Sec¬ 
ond, there is an intermediate range of energies where the 
density of states is dominated by the background value 
but Landau levels are well resolved. Third, one observes 
that the background density of states is equal to that in 
the absence of magnetic field (quadratic in energy). Fi¬ 
nally, the magnetic field creates a finite density of states 
at the degeneracy point which depends on the magnetic 
field. 


IV. CONDUCTIVITY AT CHARGE 
NEUTRALITY 


In this Section, we calculate the conductivity of a dis¬ 
ordered Weyl semimetal in the presence of a quantizing 
transversal magnetic field. Here we restrict ourselves to 
the case of weak disorder and to zero chemical potential, 
p = 0. We use the Kubo formula for the real part of the 
longitudinal conductivity. 


^ f de frie) f d^p 

a.-Au,,T)= j ^ — J 


{2tt) 


X Tr 


I G^(e,p) - G^(e,p) - uj,p)j^ 

+G^{e + uj,p)jl’^ G^{e,p) — G"^{e,p) jx'^ ■ (53) 


Here jx = evax is the bare current operator and j*’’ = 
V^'^jx is the current vertex dressed by disorder, see Ap¬ 
pendix B. The effect of disorder manifests itself in 
the replacement of bare Green’s functions by impurity- 
averaged matrix Green’s functions (8) and in the appear¬ 
ance of the current vertex corrections As discussed 
in Ref. 19, the calculation of the conductivity in Weyl 
semimetals requires taking into account vertex correc¬ 
tions even for point-like disorder (similarly to graphene). 


In the absence of magnetic field, the inclusion of ver¬ 
tex corrections away from the Weyl point leads to the 
difference between the transport and quantum (coming 
from the single-particle self-energy) scattering timesd® 
= (3/2)r*. In what follows, we first evaluate the con¬ 
ductivity without vertex corrections and then include the 
vertex corrections at the end of the calculation. 

Starting from Eq. (53), setting W’' = 1, evaluating the 
trace, and taking into account the orthogonality of wave 
functions of different LLs, we find that the Kubo formula 
for the conductivity without vertex corrections takes the 
following form in the LL representation: 


(T) = 


C de 1 eH 

J ^cosh2(^) 


dpz 

2ti 


ImGfi(e,n,p2) ImGf2(e,n,p^). (54) 


Since the self-energies for the zeroth Landau level differ 
from those for higher Landau levels, we have to distin¬ 
guish between the zeroth Landau level and the others. 
This is also true for the case of vertex corrections. For 
low temperature, T < 17, the conductivity is dominated 
by the contribution of the zeroth LL. For higher temper¬ 
atures, excitations to higher LLs are possible and there¬ 
fore, the conductivity is determined by the contributions 
of the zeroth LL, separated and overlapping LLs. 


A. Low temperatures T <C Zeroth Landau Level 

We consider first the situation when the contribution 
of the zeroth LL is dominant. This is the case under the 
following two assumptions: (i) the zeroth LL is separated 
from the first one, which is the case if the condition A < 17 
is fulfilled; (ii) the temperature satisfies T < 17, so that 
excitations to higher LLs are suppressed exponentially. 
In this case, the integral over energy e is dominated by 
the contribution of the zeroth Landau level. We note that 
for 17 A and weak disorder 7 < 7c there is no room 
for the regime A > 17, since this would imply 'yil/v^ > 1 
whereas yA/u^ < 1. Furthermore, the current vertex cor¬ 
rection for energies close to the Weyl node turn out to be 
small, U*‘'(e TC 17) ~ A/17 TC 1, see Appendix B. There¬ 
fore, in the regime of the dominant zeroth LL contribu- 
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tion, we will ignore the difference between the quantum 
scattering time and transport scattering time. 

Using ImEi ~ A and ImS 2 — 0 and disregarding the 
real part of self-energies, we get 


Gn(e,n,p^) 


G22(£,'n,Pz) 


_ " ' _ ('r.c'i 

(e + iA — vpz){e + vpz) — 

£ + iA — vpz 

{e + iA- vpz){£ + vpz) - U2(n -|- 1) ’ 

(56) 


Substituting these in Eq. (54), we arrive at 

(0) f°° de 

= 7^ / — - rTTTA Qri(e), (63) 

27r2u J_^ dTcosh^ (^) 


dz 

Qn(£) = J —lmG^^{£,n,Pz)ImG22{s,n,Pz). 


(64) 


The evaluation of the integral in Eq. (64) then yields 


Substituting Eqs. (55) and (56) into Eq. (54) and setting 
e = 0 in Green’s functions, we arrive at (z = vpz)'- 


Uxx (T) 


2tt‘^v 



(h __ 

27r [(z^ + 


X 


_ + 1 ) _ 

{[z^ + il.‘^{n+ 1)]2 -I- A'^z'^}' 


(57) 


Neglecting ^ U (the condition of separation of the 
lowest Landau level) in the denominators for higher Lan¬ 
dau levels, we see that the sum over n > 0 converges and 
gives the contribution ~ e^A^/(flv) to the conductivity, 
whereas the n = 0 term yields e^A/v. Thus, for xl ^ U 
the total conductivity at p = 0 and T <C U is domi¬ 
nated by the contribution of the zeroth Landau level and 
is given by: 


^ XX (T < U < v^/y) 


A 

(27r)2 V 


IbTT^ 


oc yiL. 


(58) 


The resulting conductivity Eq. (58) is proportional both 
to the disorder strength and to the magnetic field. 


B. High temperatures, T 3> U 


r 

Qnis) = 


v/s 2 -U 2 (n-f l)-r 2 - 2 ier 

£{2n -I- 1) -I- zT 
U 2 -h 4i£r 


(65) 


Vs^ - U2n -r^ + 2i£T 


Let us now include the vertex corrections. The to¬ 
tal conductivity is then given by Eq. (63) with the re¬ 
placement Qn —> Qn, where includes the dressing 
of the current operator by disorder (“transportization”). 
The vertex correction U*‘'(e ^ U) is calculated in Ap¬ 
pendix B: 


-I- iisT 
U 2 -h fieT' 


( 66 ) 


The inclusion of vertex corrections replaces dzeT with 
(8/3)zer in the denominator of Eq. (65). In zero mag¬ 
netic field this yields axx = {i/2t)u^xx , in agreement with 
Ref. 19. Below we will see that the effect of vertex cor¬ 
rections in magnetic field is captured by the replacement 
of t"? with T**' = (3/2 )t'^ in the Drude-like formula for 
the magnetoconductivity. 

A detailed evaluation of Qn for T ^ U is given in 
Appendix C. For e > e*, when the LL broadening is 
dominated by the background, E = 2Ae^/r2^, we have 


For higher temperatures, T ^ fl, energies e ^ are 
involved in the thermal averaging, so that we need to 
evaluate the contribution of high Landau levels to the 
conductivity. For e ^ U we neglect the difference be¬ 
tween the self-energies: Si = E 2 . As before, we include 
the real part of self-energies into the shifted energy e and 
drop the tilde everywhere. The imaginary part of the 
self-energy is written through the Landau-level broaden¬ 
ing: ImS(e) = —zr(e). 

The Green functions for £ U take the form: 


Gfi(£,n,p^) 


G§2{£,n,pz) 


£ + vpz + zT 
(e -I- zr)2 — z;2p2 _ ’ 

£ + vpz + zE 

(£ -I- zE)2 — z;2p2 _ + 1) ’ 


(59) 

(60) 


yielding with z = PzV 

T- „R „ £^ -I- z^ -I- -I- E^ -I- 2 £z 

““ (£ 2 _^ 2 _^ 2 „_r 2 ) 2 + 4 £ 2 r 2 ’ 


(61) 


ImG?o ~ -E 


£^ -f z^ -I- 0^(n -|- 1) -I- E^ -I- 2 £z 

(£2 - z2 _ f]2(„ 1) _ r2)2 4g2r2 ' 


(62) 


S*' 

n—0 


4E£4 

U2[(4£E)2 -h9U4/4]' 


(67) 


As a result, the contribution of this energy region to the 
conductivity reads: 

^ /cosh2(^) (8A£3)2+9E!8/4- 

This expression can be cast in the form of a conventional 
Drude-like formula for the magnetoconductivity (for re¬ 
cent review see Ref. 46) with the £-dependent transport 
scattering time t*’'(£) and effective cyclotron frequency 

Wc(£), 


a 


D 

XX 


Gtt 


de h'{s)T^^{e) 

4rcosh2(^) l+w2(e)[rf(£)]2- 


(69) 


Indeed, using the self-consistency relation between the 
densities of states and scattering times in magnetic field 


K£)t‘’'(£) 


3 

dzry ’ 


(70) 
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the semiclassical expression for the cyclotron frequency 
in the linear spectrum 


Wc(e) 



2e' 


(71) 


and expressing the quantum scattering rate through the 
broadening 


r(e) = 2A 




1 

2T^{e) 


3 

4T‘''(e) ’ 


(72) 


we obtain Eq. (68) from Eq. (69). 

Furthermore, as demonstrated in Appendix C, in the 
energy range 12 < e < e*, where Landau levels are sepa¬ 
rated, the conductivity is still dominated by the low-lying 
LLs with n < N^, leading to Eq. (68) also in this regime. 
It thus turns out that the semiclassical Drude formula 
(69) is valid for all temperatures T 12. A similar re¬ 
sult was obtained for graphene at the charge neutrality 
point in Refs. 47 and 48. 

Let us now analyze the magnetic-field dependence of 
the conductivity for T > 12. The denominator of the 
integrand in Eq. (68) is dominated by for e < e** = 
12(12/A)^/^. Therefore, for 12 < T < e** we obtain 

4e^ A12^ r de e® 627r^ 1 

“ 9)r2 DiF J cosh^ (#) ^ ^ W ^4124 IP' 

(73) 

For higher temperatures, T > e**, we neglect 12® in the 
denominator of the integrand in Eq. (68), which yields 

A12^ f de e® 

P vT J cosh^ (^) (8Ae3)2 27r7 ' 

(74) 

This is just the conductivity in the absence of magnetic 
field. The H dependent correction to this result is non- 
analytic in H: 

^ 12 ^/® 77 ^/® 

^ (^^8/ y/3y 

(75) 

This correction will determine the low-field magnetore¬ 
sistance. 

The summary of all regimes in the temperature- 
magnetic field plane is visualized in Fig. 6. For each of 
the regimes, the scaling of the corresponding dominant 
contribution(s) to the conductivity is shown. The con¬ 
ductivity is dominated by the Drude formula (69) down 
to the lowest Landau level, T ^ 12: 



cx 77, 

T< 12, 





1 

D < T < £** 

T» 

nf22/® 

9 2 

e 

7 


ryl/3 


( 76 ) 



FIG. 6. Behavior of the conductivity in a Weyl semimetal 
with weak white-noise disorder in dependence of tempera¬ 
ture and magnetic field. Scaling of dominant contribution to 
the conductivity in each of the parameter regimes is shown. 
Equations describing borderlines between the regimes are also 
indicated. 


In the last regime T ^ e**, the correction to the conduc¬ 
tivity is given by Eq. (75). 

As is seen in Fig. 6, the limits iJ —0 and T —>■ 0 
are not interchangeable for weak disorder. Specifically, 
setting first 77 -i- 0 and then T —0 yields 

2 2 

lim lim a^xiH.T) = -. (77) 

T^oH^o 27r7 ^ ^ 

On the other hand, if we first perform the limit T —>■ 0 
and then 77 —)■ 0, we obtain 

lim lim cra;a;(77, T) = 0. (78) 

tt-j-OT-j-O 

This demonstrates a peculiarity of the transport prop¬ 
erties at the Weyl point. As has been mentioned in the 
Introduction, a similar behavior was also found in the ab¬ 
sence of magnetic field. Ref. 13. Specifically, at 77 = 0 it 
was crucial to distinguish between the order of the limits 
T —>■ 0 and w —0 for weak disorder, while in the strong 
disorder regime, the order of the limits was interchange¬ 
able (see Refs. 12 and 13). In the presence of magnetic 
field, the limits a; —)■ 0 and T —>■ 0 become interchange¬ 
able also for weak disorder because of a finite density of 
states at the Weyl node. The reason for the elimination 
of the non-interchangeability of two limits by magnetic 
field (for any disorder strength) or by strong disorder is 
the generation of a finite density of states at the Dirac 
point (e = 0). 


C. Conductivity for strong disorder 

We now briefly discuss the magnetoconductivity in the 
case of strong disorder, /3 = 'j/jc ~ yA/n® 1. As dis¬ 
cussed in Sec. IIIC, in this case the broadening of all 
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Landau levels is of the order of the ultraviolet cut-off: all 
Landau levels overlap. To calculate the conductivity, we 
use the semiclassical Drude formula (69) complemented 
by Eqs. (46), (48), (70), and (71). The parameter gov¬ 
erning the magnetic-field dependence of the conductivity 
now reads 

o 2 o 2 

cj.ieHe) (79) 

Here we have replaced the energy e by the level broaden¬ 
ing r ~ A. The main contribution to the conductivity is 
independent of magnetic field, whereas a weak magnetic 
field yields a quadratic-in-correction (here we do not 
write numerical prefactor in the i7-dependent correction) 

I" de 1 

Ti + nyA^^^[ A^J’ 

g2o4 

-oc -i7^. (80) 

For strong disorder, the i7-dependent conductivity cor¬ 
rection vanishes quadratically. This should be contrasted 

with the conductivity for weak disorder which is non- 
analytic (proportional to 


V. CHARGED IMPURITIES 


A. Screening 


In the preceding part of the paper we considered a 
model of white-noise disorder. We are now going to gen¬ 
eralize the obtained results onto the more realistic case 
of screened Coulomb impurities. The potential of such 
an impurity is given by 


C/(k) 


47re^ 

eoo(fc^ + ^ 2 ) ’ 


(81) 


where e^o is the background dielectric constant. The pa¬ 
rameter K is the inverse Debye screening radius. 


in order to find the parametric dependence of the con¬ 
ductivity (without numerical prefactors). 

The correlator (84) corresponds to a white-noise disor¬ 
der whose strength depends on magnetic field and tem¬ 
perature, 

r, 

T, 

(85) 

where Aimp is the density of impurities. In the limit 
T,H —>■ 0, Eq. (85) yields a divergent disorder strength 
which implies the necessity of a self-consistent treatment 
of the impurity screening. Specifically, when at 

max(fl, T) ~ Eimp = AVp'^ ( 86 ) 

the quasiparticle broadening 7 ( 77 , r)max(T^, Vl?)/v'^ be¬ 
comes of the order of max(n,T), the screening will be 
determined by the impurity-induced density of states, 
yielding 

7(77, T) ~ 70 = TVj-f « 7 ,. (87) 

The weak-disorder approach is applicable under the con¬ 
dition max(f2,T) > eimp- 

Below we employ the results of the previous sections 
to the Coulomb case by replacing 7 with 7 ( 77 , T). Using 
Eq. (85), we express the condition for separation of the 
zeroth LL, A < U, as 

0 > u^/70 = eimp- (88) 

We see that for screened Coulomb impurities the increase 
of magnetic field favors the Landau quantization, as op¬ 
posed to the case of white-noise disorder. In particular, 
in the limit 77 —)■ 00 , the zeroth LL is always separated. 
This demonstrates a crucial role played by the 77 depen¬ 
dence of the screening. 


7(77,T)=^imp(|; 





B. Conductivity 


2 47 re^ dn f ^ T, 

Coo TTSooV^ I tt^T^/S, U < T. 

Here dn/dfi is the fermion compressibility and we have 
neglected the effect of disorder on the thermodynamic 
density of states. In the limit T, 77 —)■ 0 one should in¬ 
clude the impurity contribution in a self-consistent way, 
as we will discuss below. 

In the following, we shall assume that the “fine- 
structure” constant is not small, 

e^/v > 1. (83) 


In order to find the conductivity, we should substitute 
7 ( 77 , T) ^ 7 (C^[max(U, r)]“^ for 7 in the results for the 
conductivity of a system with white-noise disorder. Since 
we do not keep numerical factors, we also disregard the 
vertex corrections (that only modify these factors, see 
above). We begin by considering the low-temperature 
regime, T ^ U, when the screening is controlled by the 
magnetic field. If U < £imp, all LLs overlap, and the 
conductivity is essentially equal to that at zero magnetic 
field, see below. In the opposite situation, U > eimp, 
only the zeroth LL contribute. With A ~ AimpU^/U^, 
this yields 


In a realistic situation it is of order of unity, so that k is of 
the order of characteristic values fctypicai ~ max(f 2 ,r)/u 
of the wave vector k. Condition (83) allows one to 
describe the screened Coulomb disorder by an effective 
pointlike correlator 

{U{r)U{r'))^j{H,T)6{r-v') (84) 


^XX 


e^A 

V 


^ ‘^imp 


e^AimpU^ 1 

77’ 


(89) 


This result agrees with the result obtained by Abrikosov 
in Ref. 39. We observe that the dependence of conductiv¬ 
ity on the magnetic field for Coulomb impurities differs 
strongly from that for white-noise impurities. 
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FIG. 7. Behavior of the conductivity in a Weyl semimetal with Coulomb impurities in dependence of temperature and 
magnetic field. Scaling of dominant contribution to the conductivity in each of the parameter regimes is shown (the Fermi 
velocity v is set to unity). Equations describing borderlines between the regimes are also indicated. The symbol ^ denotes 
numerical coefficients. 


We turn now to the domain T > characterized by a 
temperature-dependent screening. Here contributions of 
higher LLs are important. We use the semiclassical ex¬ 
pression (69) for the conductivity together with Eq. (71) 
and 


1 ^ l{T)e^ ^ 

T(e) 


(90) 


A closer look reveals a necessity to distinguish between 
several regimes. For e < one can neglect 

the unity in the denominator of the integrand in Eq. (69). 
Therefore, for T > the conductivity of a system 

with screened Coulomb impurities takes the form 


For weaker magnetic fields, fl < e?„„p/r^, the energies 
dominating the magnetic-field dependence of the conduc¬ 
tivity are below eimp. As a consequence, one should re¬ 
place £ by Einip in the parameter WcT: 

WcT ~ • (94) 

^imp 

The iif-dependent correction to the conductivity for < 

is 


g2 047^11 

Scr^x ~-14— OC -H'^. (95) 

^^imp 


vVL^ ^ ij2 • 


(91) 


For T < the conductivity is dominated by the 

range of e where we can neglect the magnetic field: 


e^T'^ 

7(0 ■ 


(92) 


The i7-dependent correction to this result turns out to 
be different in the two subregimes. The first subregime 
is defined by < VL < The result is 

obtained similarly to Eq. (75): 




e 02 112/37.1/3 
7(T) £imp 


527.13/3 

imp 


112/3 


cx 


( 93 ) 


We thus see that the non-analytic (77^/®) magnetoresis¬ 
tance does not survive the limit iJ —>■ 0 for the case of 
charged impurities, in contrast to the case of white-noise 
disorder. 

For the lowest magnetic fields H < eimp and tempera¬ 
tures T < Einip, all Landau levels overlap. We thus ex¬ 
pect that a zero-77 calculation, as recently carried out in 
Refs. 26 and 27, should be applicable. In our notations, 
the result of those papers is written as 


-N, 


imp 


‘-'imp 

V 


(96) 


In this regime, disorder is strong. This means that the 
77-dependent correction to the conductivity can be cal¬ 
culated under the assumption that l/r ~ T .... Eimp and 
ujc ^ n^/r. The 77-dependent correction to Eq. (96) for 
H r <C Eimp is calculated similarly to Eq. (80) and 
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takes the form: 


d(Tx 




cx -i?" 




(97) 


Finally, in the regime T ^ £imp an analogous con¬ 
sideration with l/r ^ n^/fTimp yields no H dependent 
correction to the leading order. 

The conductivity in all regimes is summarized in Fig. 7. 
For each regime, the scaling of the dominate contribu¬ 
tions is shown. We can summarize the results according 
to the parameter T/eimp, yielding to 


^ ^imp 


€ SiT] 


1 


1 - 







(98) 


for T/eiinp < 1. In the opposite limit, we need to distin¬ 
guish between more regimes 


^ ^imp 


1 

T2 


e^e? 

^imp 


g2j.4 / 


^£fmp 


1 - 


1 - 


^^2/32^1/3 


-imp 


fl>T, 

•ill 

3 3/2 

fiUlP ^ O ^ 


ve'. 


imp 


j 


-11 

•"imp 




"■imp 
2^2 ■ 


(99) 

Here and in Fig. 7, the symbol ^ denotes numerical co¬ 
efficients. We observe that the limits T —>■ 0 and H —>■ 0 
are interchangeable for Coulomb impurities. 


VI. MAGNETORESISTANCE 

Using the results of the previous sections, it is straight¬ 
forward to evaluate the magnetoresistance described by 

Ap{H) = ( 100 ) 

Pxx\y) 

The resistivity is given by pxx = o'xxiio'xx + ^xy)- Since 
we considered the case of zero chemical potential p = 0, 
the Hall conductivity is zero and we can calculate the 
magnetoresistance with 

Ap(i7) = - 1. (101) 

^XX\^ ) 

In the following, we calculate this relative magnetoresis¬ 
tance for both models of disorder, pointlike impurities 
and charged impurities. 


A. Pointlike Impurities 

We fix the value of T 7 and analyze the evolution of the 
magnetoresistance with increasing magnetic field (which 



FIG. 8. Dependence of the magnetoresistance on magnetic 
field (D^ oc H) for the case of white-noise disorder (weak dis¬ 
order regime). Magnetoresistance in the range of relatively 
weak magnetic fields, Q, <C is shown. For larger fields 

the magnetoresistivity would vanish. Scaling of the magne¬ 
toresistance in different regimes, regimes boundaries and the 
corresponding values of Ap are indicated. 


corresponds to a vertical cross-section in Fig. 6 ). Using 
Eq. (76), we obtain the magnetoresistance 


Ap = < 


7r4/3 uU2/3 

122/3 ^l/3y’ 

189 ^ 

1247r4 72r6 “ ’ 


2I/22-3/2 


^ i;3/2 

2l/2y3/2 


„3/2 


< U < T, 


Stt^ - 1, T <C U <C—. 

7^SU j 

( 102 ) 

For U > ^ 3 / 7 , the magnetoresistance vanishes, because a 
further increase of magnetic field leads again to a regime 
of overlapping LL. The magnetoresistance in different 
regimes is visualized in Fig. 8. 

It is worth emphasizing that the vanishing of the den¬ 
sity of states of Weyl semimetal at e = 0 in zero H 
translates into the non-analytic, behavior of the 

magnetoresistance for weak pointlike impurities, which 
persists down to H = 0. This should be contrasted 
with the i73/2_jjjg^gjjgiQj.ggjgig^jj(.g graphene found in 
Ref. 47. The square-root magnetoresistance in graphene 
does not actually survive the limit 77 —>■ 0, since even 
weak white-noise scalar disorder is marginally relevant 
in graphene and establishes a finite density of states at 
the Dirac point. As a result, the true 77 —>■ 0 asymp¬ 
totics of the magnetoresistance in graphene is parabolic 
(although the crossover between 77^ and 77^/^ may occur 
at very weak magnetic fields). By contrast, weak white- 
noise disorder in 3D systems with linear spectrum is not 
capable of establishing a finite DOS at £ = 0. This leads 
to the remarkable observation of the non-analyticity of 
the transversal magnetoresistance. This non-analyticity 
can be considered as a finite-77 counterpart of the non¬ 
commutativity of the limits T —>■ 0 and a; —>■ 0 that takes 
place in zero magnetic field. 

We also stress that in the limit of low temperatures 
T ^ 0, the magnetoresistance for pointlike impurities is 
very sharp, with the maximum at 77 oc —>■ 0 growing 
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FIG. 9. Magnetoresistivity for a Weyl semimetal with 
Coulomb impurities in the limit of low temperatures, T < 

S^imp • 

as Ap cx 1/T^ — oo. This is a manifestation of the non¬ 
commutativity of the limits T —0 and i/ —0 discussed 
at the end of Sec. IV B. 


B. Charged Impurities 


For the case of charged impurities, we perform the 
analysis in a similar way. Specifically, we fix the pa¬ 
rameter T /Eimp and then consider the magnetoresistivity 
when the magnetic field is swept. 

For T/^imp < 1, Eq. (98) is used to obtain the magne¬ 
toresistivity, reading 


Ap 




^imp 



^imp 


^ ^imp: 

£imp ‘C fl. 


(103) 


This type of behavior in the high field limit was identi¬ 
fied in an early work by Abrikosov, Ref. 39. For lower 
magnetic fields, the LLs overlap, which implies that the 
conductivity is essentially the same as in the absence of 
magnetic field. Therefore, the magnetoresistance van¬ 
ishes quadratically, as visualized in Fig. 9. 

In the temperature limit T/e-irap > 1, we use Eq. (99) 
to obtain the magnetoresistivity. This leads to 


A„ 


imp 

^2/3tI/3 


-imp 


imp 


- 1 , 


- 1 , 


"imp 




y2 


3 3/2 

fiinp ^ p ^ 

3/2 


(104) 


In the high-field limit, the magnetoresistivity is again 
linear in magnetic field. The low-field limit is dominated 
by a quadratic dependence on magnetic field. This is 
visualized in Fig. 10. 



FIG. 10. Magnetoresistance of a Weyl semimetal with 
Coulomb impurities for sufficiently high temperatures, T > 
Simp, in the whole range of magnetic fields (a). The low-field 
region is magnified in panel (b). 


C. Comparison to experiment 

As has been mentioned in the Introduction, two re¬ 
cent works^®’^^ reported a strong, approximately linear 
magnetoresistance in the Dirac semimetal Cd 3 As 2 in the 
range of strong magnetic fields. We have found that a 
model of pointlike impurities yields a very different be¬ 
havior of the magnetoresistance and thus cannot explain 
the experiment. On the other hand, a more realistic 
model of disorder—Coulomb impurities—does lead to a 
linear magnetoresistance in the region of strong magnetic 
fields, yielding a tentative explanation of spectacular ex¬ 
perimental observations. In particular, the experiment 
of Ref. 36 shows a linear longitudinal resistivity starting 
at the magnetic field = l.IT for T = 300K. This is 
in a good agreement with the onset of linear magnetore¬ 
sistance at D > T [Eqs. (103) and (104)] which corre¬ 
sponds to Ff ^ IT for T = 300K. At the same time, in 
the range of weaker magnetic field, H < VT the experi¬ 
mentally measured magnetoresistance is parabolic, which 
again agrees with our result for Coulomb impurities. The 
slope O.dmDcm of the linear resistivity in experiment^® 
agrees with our estimate ~ O.OmDcm obtained under the 
assumption that the carrier concentration is of the or¬ 
der of density of Coulomb impurities. Furthermore, the 
magnetoresistance in Ref. 32 shows two inflection points 
similar to our result of the magnetoresistance visualized 
in Fig. 10 (b). 
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Let us point out that we have performed the analysis 
for the neutrality point, /i = 0, while the experiment is 
carried out at a finite carrier density. The latter leads to 
an observation of Shubnikov-de Haas oscillations which 
are absent in our /i = 0 theory. We expect that our 
theory should become applicable to a system with finite 
carrier density for sufficiently strong magnetic field, such 
that the chemical potential is located between the zeroth 
and the first LL. 


VII. SUMMARY AND DISCUSSION 

To summarize, we have developed a theory of the mag¬ 
netoresistivity of Weyl and Dirac semimetals at the neu¬ 
trality point. We have considered two alternative models 
of disorder: (i) short-range impurities and (ii) charged 
(Coulomb) impurities and treated the impurity in the 
framework of the self-consistent Born approximation. We 
have found that an unusual broadening of Landau levels 
leads to a variety of regimes of the resistivity scaling in 
the temperature-magnetic field plane. 

The behavior of the magnetoresit ance is essentially dif¬ 
ferent for two types of disorder, with the difference orig¬ 
inating from the dependence of screening of charged im¬ 
purities on magnetic field. In the limits of strongest mag¬ 
netic fields H, the magnetoresistivity vanishes as 1/H for 
pointlike impurities, Eq. (102), while it is linear and pos¬ 
itive in the model with Coulomb impurities, Eqs. (103) 
and (104), in agreement with experimental observations 
in Refs. 36 and 37. The prefactor of this linear magne¬ 
toresistance is approximately temperature-independent 
for low temperatures. In the low-field limit, we find a 
quadratic magnetoresistivity for screened Coulomb im¬ 
purities, Eq. (103). By contrast, the low-field magne¬ 
toresistance for pointlike impurities shows a non-analytic 
behavior, oc 

It should be emphasized, however, that our theory was 
developed for the Dirac point, /i = 0, whereas the ex¬ 
periments have been carried out at a finite carrier den¬ 
sity which manifested itself in a non-zero Hall resistivity 
and in Shubnikov-de-Haas oscillations of the resistivity. 
For non-zero fx the present theory is justified only in the 
range of sufficiently high temperatures or magnetic fields, 
^ < max[f2,T], for which the chemical potential is lo¬ 
cated between the zeroth and first Landau levels. 

A generalization of the theory on the case of finite /x 
and lower magnetic fields remains a prospect for future 
research. Another interesting generalization of our the¬ 
ory would be the analysis of the effect of magnetic field 
on the critical region around the transition (7 = 7c) be¬ 
tween weak and strong disorder. 

Further, our theory assumed the absence of internodal 
scattering between different Weyl points and diagonal 
scalar impurity potential. The internodal scattering is 
believed to be weaker than the intranodal one in real¬ 
istic Weyl semimetals. However, even weak internodal 
scattering is crucially important for quantum effects, es¬ 
pecially for those related to the famous chiral anomaly. 


Since the transverse magnetoresistance (in contrast to 
the longitudinal magnetoresistance) is not directly af¬ 
fected by the chiral anomaly, we expect that our semi- 
classical calculation based on the SCBA will not change 
qualitatively in the presence of internodal scattering, pro¬ 
vided the two Weyl points correspond to the same energy. 
In the case of non-degenerate Weyl points, the intern¬ 
odal scattering is expected to establish a finite density of 
states at any energy even in the limit of zero magnetic 
field. In this situation, we expect that the non-analytic, 
behaviour of the magnetoresistivity at iL —> 0 for 
weak pointlike impurities (stemming from the vanishing 
DoS) would be smeared at the lowest H. A detailed 
study of this case is a prospect for future work. In a 
similar way, the nondiagonal matrix elements of disorder 
potential should not affect our SCBA results qualitatively 
as long as they do not induce a finite DoS. 

Finally, the Coulomb interaction between quasiparti¬ 
cles and the particle-hole recombination in a finite geom¬ 
etry are additional sources of the magnetoresistance at 
the charge neutrality point, similarly to 2D compensated 
systems, see Refs. 49, 50 and references therein. These 
effects are expected to be important for the description 
of the high-temperature magnetoresistance. 

Note added. After the completion of our paper, we 
became aware of very recent related works on magneto¬ 
transport in Weyl semimetals. In Ref. 51 the magnetore¬ 
sistance is analyzed for a model with smooth random 
potential. In Ref. 52, a transversal magnetoresistance 
in Weyl semimetals with Coulomb impurities is calcu¬ 
lated for a small “fine-structure constant” [cf. Eq. (83)], 
and the result is consistent with our results obtained for 
^ jv > 1. In Ref. 53, the role of Fermi arcs in a finite- 
width slab of a Weyl semimetal subject to a static mag¬ 
netic field is addressed in the context of non-local re¬ 
sponse. While our calculations have been performed for 
macroscopic samples, it would be interesting to investi¬ 
gate a possible effect of Fermi arcs on the transversal 
magnetoresistance in a thin-slab geometry. 
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Appendix A: Shape of the broadening of separated 
Landau levels 

In this appendix we analyze the shape of the LL 
broadening at e < e* = 0(U/A)^/^ (when LL are well 
separated). The maximum of r(£) around £ ~ Wn 
is located at £ ~ and is given by 
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order of A{e/^Y {zeio-H result). 

On the left side of the peak, for e < Wn — Fat, the 
solution of the self-consistency equation (40) yields 


r(e) ~ 2A 


02 


A\/Wjy 

[2{Wn - e)]^/- AvWw' 


(A2) 


which matches Fat at e ^ Wn — Fat — Ae^ j'A}. The 
decrease of left side of the peak from F(e) ~ Fat to F(£) ~ 
Ae^ jA is thus very sharp. The Landau-level broadening 
for the case of well separated levels is shown schematically 
in Fig. 11. 


FIG. 11. Schematic illustration of the Landau-level broad¬ 
ening. At energies e < £i, such that |£i — Wa^I ~ Fiv, 
the broadening is given by Eq. (A2). The maximum of 
ImE is Ftop — FAr3^^^/22/^ ~ I.IFat and is achieved at 
£2 = Wn -|- F n 12^^^ ~ VFjv -I- O.GSFiv. At £3 ~ VFjv -I- Fat the 
tail (Al) develops. At £4 ~ Wn + 0(0/£)^ this tail reaches 
the background 2A[N -\- 1). 


Appendix B: Vertex corrections 

In this appendix, we analyze the vertex corrections in 
the diagram for the conductivity. We calculate the fol¬ 
lowing diagram that describes the dressing of a current 
vertex by disorder lines: 


Ftop = F(W^)3i/2/22/3, where F(W^) = {A/2f/^wU\ 
For brevity, we use the abbreviation Fat = F(lFAr). 

For £ > VFat -f Fat, the peak in F(£) decreases as 


F(£) ~ Fat 


2F 


N 


e-Wn 


A I Wn 

71 V e-Wn' 


(Al) 


and reaches the value of the background at £ ~ Wat -I- 
fl(n/£)^. This value is always smaller than Wn+i — 
Wn + /{2e). Thus, in the range Wn + fl(n/£)^ < £ < 

LFat+i — F(WAr+i), the Landau level broadening is of the 



X G^{E,Py,P:,,X,x')aa;G^{E,Py,p^,x',x). 

(Bl) 

The diagonal elements for the resulting matrix vanish 
after the integration over x, x' and py because of the 
orthogonality of the Hermite functions. The off-diagonal 
elements Vx,i 2 are equal. They read 


r ^ jeH _ £ - Sf -f vp^ _ £ - - vp^ _ 

J 2tt 2'kc {e— Y.f'— vp^){e— + vp^) + A^n {e— — vp^){e— + vpz)+ A^{n + \) 


For energies away from the Weyl point, the difference between the two self-energies, Ei and S 2 , can be neglected. 
Furthermore, for high Landau levels we can express the Landau level index n in terms of the momenta in px and py 
direction. The correction to the current vertex then simplifies as follows 


Vx,i2 =en7 


dp^ {e — Y^){e — Y^) — v'^pI 

(27r)3 ((£ - E^)2 - n2p2)((e _ e«)2 - ^ 2^2 


fl2)’ 


(B3) 


which can be cast in the form: 

Vx,12 =- 


(£ — S"^)(£ — E-^) enq 

f 

1 

1 

(£ - S«)2 - (£ - S'4)2 + 02 TT j 

27rP 

_(£-E^)2-n2p2_^|2 

{e — S-^)2 — n2p2 


+ 


eny 

Btt 


dp, 
27r^ 


1 1 

{e — S'^)2 — n2p2 _ J 72 ("g _ _ ^ 2^2 


(£ — Y"^)^ + (s ~ Y^)^ — A^ eny 
"(e- E«)2 - (£ - Y^y +~A^^ 


' dp 2 
2^^ 


1 1 
{e — E'^)2 — t|2p2 — J 72 _ 5]fl)2 _ y2p2 


(B4) 


The integrals over momenta here can be identified with those for the self-energy. The shift in the denominator by 
A of the integrands is unimportant (similarly to the difference between Si and E 2 ) and can be neglected; it is then 
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sufficient to keep in the prefactors of the integrals. The corrections to the current in cc-direction then reads 


/ (£-E41)(£-E«) 

■ 

Y^ ' 

1 

6 

r E^i E^ 1 

1 

V(£-E«)2- (e-E-4)2_^f^2 
1{£-Y^Y + {£-Y^Y 

_£-E41 

r E^ 

£-YR_ 
Y^ ■ 

_£-E^ £-E«_ 

6 (£ - E«)2 - (£ - E4l)2 + f 72 

'cn 

1 

M 

£-YR_ 

)■ 



For weak disorder, using <C e, after some algebra we arrive at: 


(B5) 


ev AieT 

^ y 4ier + 112 ■ 


(B6) 


This result matches to the result in Ref. 19 for the case H = 0. The summation over the ladder of impurity lines 
yields 




1 AieT{e) 

3 4i£r(e) + Q? 


(B7) 


leading to Eq. (66) and hence 

y (s) = 


(B8) 


in the Drude formula (69). 

This result is not applicable at the Weyl point, where Ei E 2 . To simplify Eq. (B2), we can use E 2 = 0 which is 
justified for weak disorder. Equation (B2) simplifies then as follows 


f dpz jeH e + vpz 

J 2tt 27rc {e — — vpz){e + vpz) + {e — 


e-Y^- vpz 

vpz){e + vpz) + fl^(n + 1) 


(B9) 


We find that this integral is proportional to A/Q, and hence the vertex corrections are small. This is in agreement 
with the neglect of vertex corrections in the calculation of the strong-77 conductivity dominated by the lowest Landau 
level in Ref. 39. 


Appendix C: Details of calculation of the conductivity for high temperatures 


In this Appendix, we present details of the evaluation of the conductivity in the regime of high temperatures when 
many Landau levels contribute to the result. We first calculate the conductivity without current-vertex corrections 
and include the vertex corrections in the end of the calculation. Starting from Eqs. (63) and (64), we represent Qn{^) 
with = 1 as a sum of the two terms related to ImG(\ and ImG |2 Eq. (64), respectively: 

dz +Vt^n + T‘^ A-2 £z){£^ + + Vi^{n+l)+ T‘^+ 2 £z) - Ae^z'^ _ p(i) (odi) 

y [(£2 _ ^2 - - r2)2 -H 4£2r2][(£2 _ ^2 _ + 1) _ r2)2 + 4£2r2] “ , i j 

where (with = fl^n + T^) 


Qn{e) = r^ J 


= r 

= r 


2r2[fi2 _ 4^2(2n -f 1)] -f [4r2 -t fi2(2„ +1)] y -w^ + 


2[(4£r)2 -f _ 44 ^ 2)2 + 4e2r2^£2 -W^ + V^£2 - W2)2 A£^ 

_2r2[L!2 _ 4^2 (2„ + 1)] + [4r2 + 1^2 (2„ + l)] ^^2 _ ^2^^ ^ ^(£2- wyi)2+4£2r2 


(C2) 


2[(4£r)2 + n4]^(£2 _ vp2^^)2 + 4£2r2W£2 - ^{£‘^ - + A£-^T^ 


(C3) 

Since the term is dominated by £ ~ lF„+i, it is convenient to shift the summation over n for this term. For 
definiteness, we use the hard cut-off for the summation over Landau levels, such that the highest Landau level involved 
in the summation is (A^max — 1) + 1 = IVmax from 


Afmax-l 


N —1 

-‘’’max 




n—O 


n—1 


^ JV —1 ’ 

-‘''max -1- 


(C4) 
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where 


and we write 


n® = 

2[(4£r)2+ll4]’ 

(II) 

^iV„ax-l ll[(4£r)2 +04]’ 
n(i) + =_ 'J'^Te _ 

+ V„-l 

Qn — Qn,l Qn,2- 


Here we have split into two parts related to the asymmetry of each Landau level: 


^n,i — n 


-W^ + y(e2 - ^^2)2 + 4£2r2 

7P^WF+4pp ' 


(C5) 

(C 6 ) 


(C7) 

(C 8 ) 


(C9) 


qn,2 = 



+ wl-^ 


1 ^2)2 +4^2r2 


,Jw^-e^ + y (£2 - 11^2)2 + 4£2r2 

y(£2 - W'2)2 + 4£2r2 


(CIO) 


For weak disorder, r(£) <C £, and :+ ( 7 „_ 2 - We note 
that these functions have resonant structure and take 
their maximal values when 


£^ ~ O^n + 



(Cll) 


Comparing Eqs. (C9) and (CIO) at resonances, we see 
that the maximal value of is dominated by qn,i'- 


3/2 

gn.l ^ ^^pl /2 ’ g «,2 y2rl/2£l/2 < g„4. (C12) 

One can check that i ;+ g„ 2 in the whole range of 
energies Wn < e < W„+i. Therefore, in what follows we 
will disregard the contribution of qn, 2 , using g„ ~ gn,i- 
Comparing Eq. (C9) and (33), we find [see also Eq. (40)]: 


qn — ^ Tl 


V 2 r(")(£) 


£l£ 




(C13) 


Let us now consider energy e located between the Lan¬ 
dau levels N and fV + 1. We split the sum over Landau 
levels in Eq. (C4) as follows: 


— 1 




v^r£ 


(4£r)2 + 04 


y 2 r£ 


' N-1 


A^max -1 

E 

n—l 

A^max -1 


(4£r)2 + 04 


E! Qn + qN + qN+1 + E! Qn • 
n=l n=N+2 J 

(C14) 


When the Landau-level broadening is smaller than the 
distance between the neighboring levels, Wm — Wn+i ~ 
0 ^/£ > r(£), i.e., £<£** = 0(0/x4)^/^, we can neglect 


r in the contributions of all Landau levels with n < N. 
Replacing the sums by integrals, we find 


AT-l 

E 

n—l 


'/2n^n 

■— \/£2 — 02 n 

n—l ^ 


4v^ £3 

“^ 02 ’ 


(C15) 


This contribution to Eq. (C4) is by a factor £^/0'^ larger 
than the n = 0 term, Eq. (C5), so the latter can be 
neglected. For n > iV + 1 we expand q^ in F and get 


E 

n^N+2 


4''max J- 

E 

n^N+2 


v^r( 2 £^ - O^n) 

( 02 n — ^2^3/2 


2 v^£r 




2v^£3r 


02y02(lV + 2)-£2 


. (C16) 


The term proportional to -s/lVmax exactly cancels the 
contribution of Eq. (C 6 ) in the sum over n, Eq. (C4). 
The second term here is of the order of F/O [remember 
that O^iV < £^ < 0 ^(A^ + 1 )[ and can be neglected for 
£ < 0(0/x4)3/2 compared to the contribution oi n < N. 
Then Eq. (C4) takes the form 


Ec. 

n =0 


2r£4 

02[(4£r)2 +^4] 


■4 ^^2 

.3 + 71^^'^^ + '^^+'^, 


(C17) 

According to Eq. (C13), the second term in Eq. (C17) is 
given by 


02 

7!7 


(giv + gAT-l-l) 


02 /r(^ 
£ \ A£ 



(CIS) 


This term may dominate the sum in Eq. (C17) only for 
£ < £* ~ 0(0/A)^/3 and when £ is located close to a 
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Landau level center. We will return to this case later 
and first analyze the opposite (simpler) regime of £ > e*. 

For £>£** = the level broadening in 

is larger than the distance \e — Wnj for Landau levels 
sufficiently close to N: 

Fjv > |£-W„| ~ \Wn-W^\ 

For such values of n we can neglect as compared 

to eF in Eq. (C9): 


Qn — 


V^' 


(C20) 


In this case, the sum over n in Eq. (C14) is written as 


Afma,x-1 (N-Nr-1 N+Nr 
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02 ■ 


(C21) 

(C22) 

(C23) 


We see that for £ > £** again the first term {n < N — 
Nr) dominates, yielding the same result as for £ < £** 
(clearly, iVp <C N, in view of r(£) ^ £). Including the 
vertex correction calculated in Appendix B, we arrive at 
Eqs. (67) and (68) of the main text. 

Let us now return to the case of lower temperatures, 
O < T < £* ~ 0(0/A)i/3. In this case, the contri¬ 
bution of the Landau level N closest to the energy £ 
in the sum in Eq. (C17) should be analyzed. In or¬ 
der to estimate this contribution, we replace the inte¬ 
gral over energy in Eq. (63) by a sum over regions of 
width Fjv around Landau levels, use Eq. (C13), replace £ 
by Wn, and replace r’^^)(£) there by its maximal value 
r(^)(Wv) = Fat ~ A2/3fIi/3]vi/6, As a result, we get 


(T 


(N) 

XX 


ATv -E ^ 

Af<(T/n)2 

WN ~ir 


r 2 w 2 

{AWnTnY + 


(C24) 


(C25) 


In addition to this contribution, there is a contribution 
of the tail at Fat < e — Wn < see Fig. 11 

and Eq. (Al). The integral over |£ — WatI is logarithmic 
since r^(£) decays as (£ — IF^v)”^ in this range and thus 
enhances the result (C25) by a logarithmic factor. We see 
that the contribution to the conductivity of the Landau 
level N is smaller than the semiclassical contribution, 
Eq. (73), by factor N/N <C 1 (up to the logarithm) and 
can be neglected. 
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